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CN ■ We construct, using simple geometrical arguments, a Wigner function denned on 

a discrete phase space of arbitrary integer Hilbert-space dimension that is free of 
; redundancies. "Ghost images" plaguing other Wigner functions for discrete phase 

spaces are thus revealed as artifacts. It allows to devise a corresponding phase-space 
^ . propagator in an unambiguous manner. We apply our definitions to eigenstates and 

propagator of the quantum baker map. Scars on unstable periodic points of the 
corresponding classical map become visible with unprecedented resolution. 
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1 Introduction 



Since the conception of the Wigner function in 1932 [1], phase-space repre- 
sentations of quantum mechanics have found a wide range of applications. 
Notably in complex quantum dynamics, they have become popular as they 
allow to compare classical and quantum dynamics on the same footing. The 
Wigner function in particular has become the standard tool whenever the full 
information contained in the density matrix, including coherences, has to be 
represented. 

In many cases the "quantum phase space" underlying the Wigner function 
is fundamentally different from the corresponding classical one in that it is 
discrete. This occurs, for example, if both classical phase-space variables are 
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cyclic, i.e., if the classical phase space has the topology of a torus so that the 
periodicity of either variable implies quantization of the other. Likewise, it 
is inevitable once, e.g., in numerical implementations of wave-packet scatter- 
ing, the spatial variable is reduced to a discrete lattice and at the same time 
restricted to some finite total length, with periodic boundary conditions. Fi- 
nally, discrete phase spaces allow to define a perfectly self-consistent quantum 
mechanics living on a finite-dimensional Hilbert space [2]. 

A straightforward generalization of the original definition to the discrete case 
leads to a Wigner function that contains the information of the density oper- 
ator in several replicas. This redundancy is reflected in so-called ghost images 
of features standing out in the quantum state [3,4]. Modified definitions that 
avoid these redundancies have been proposed in [5,6]. We here start from the 
version in Ref. [6], giving a transparent derivation in simple geometric terms 
that readily generalizes to arbitrary (not necessarily odd or even prime) integer 
Hilbert-space dimensions. It allows to construct a propagator on the discrete 
phase space, which can be defined unambiguously only for a Wigner function 
free of redundancy. As an illustrative example, we apply it to the quantum 
baker map [7]. A more detailed account of our findings will be published else- 
where [8]. 



2 Non-redundant Wigner function and propagator 

The original definition of the Wigner function [1], W(p, q) = h~ l dq' e~ ipq '^ h 
p(l + q'/2, q — q' /2) for a density operator p consists of two steps, a transfor- 
mation to sum-and-difference coordinates q — (gi + g 2 )/2, q' = q 1 — q 2 of the 
density matrix p(qi,q2) = (Qi\p\Q2), followed by a Fourier transform along q' . 
The difficulty in applying this to a discrete spatial basis \n), n — 0, . . . , N — 1 
arises from the fact that in (q, g')-space the original square lattice appears as 
"base-centered cubic" , with even and odd points in terms of n\ ±ri2- Augment- 
ing it to become "simple cubic" again by introducing fictitious extra lattice 
points where p(n,n') = lead, upon Fourier transformation (of size 2N), to a 
fourfold redundancy of the resultant Wigner function [3,4], call it V4 / doubie( m ) ) 
abbreviating the discrete phase-space vector m = (A, n). We remove the re- 
dundancy as follows: By an additional Fourier transform (size 2N), now in the 
g-direction, the redundancy takes the explicit form of an identical repetition, 
up to sign changes, of the coefficients in four N x iV-blocks within the total 
2N x 2N "Brillouin zone" . We cut out a single N x iV-block, chosen centered 
around the origin so as to preserve twofold symmetries that reflect the reality 
of the Wigner function, by multiplying with a two-dimensional discrete box 
function. Upon undoing the additional second Fourier transform (now of size 
N only), the Wigner function is recovered, free of redundancies. 
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Cutting out the N x iV-block from the Fourier-transformed Wigner function 
^double is equivalent to a convolution of W^oubie with a Fourier-transformed 
box function S 



N-l 



W(m)= Yl W douhlc (2 m -m')S(m'), 

\'n'=-N 



S(m) = s x s n , s k = < 



hmod2N k even, 

sin(7rfc/2) ^—vKk/2N fc qM 
, N sm(nk/2N) 



(1) 

(2) 



Apart from revealing "ghost images" as artifacts and removing them, this pro- 
cedure becomes indispensable once a propagator of the Wigner function is to 
be defined in an unambiguous manner. The reason is that with a redundant 
Wigner function, there is no unique way of relating the coefficients of the final 
Wigner function to those of the initial one. Denoting the transformation out- 
lined above symbolically as Tw, so that W(m) = T\yp(n), with n = (ni,n 2 ), 
we find for the discrete Wigner function propagator 

K w (m", t"; m', *') = T w K(n", t"; n', t')T^\ (3) 



where K(n", t"; n', t') propagates the density matrix from time t' to t", i.e., 
p(n",t") = E' n K(n",t";n',t')p(n',t'). 

If the time evolution is unitary, generated by a Hamiltonian H, the prop- 
agator can be expressed in terms of energy eigenstates H\a) = E a \a) as 

K(n",t";n',t f ) = Y, a ,p ^ i{Ea - E0){t "~ t ' )/h PaA Q '')Po'A' a ')^ denoting p Q ^(n) = 
(a\ni) (n 2 \P). Defining Wigner eigenfunctions as ^^(m) = T w p aij a(n), we 
have for the Wigner propagator 

K w (m", t"; m', t>) = N' 1 ^ B-^- E ^-^ h W^{m^W^{m'). (4) 

a,f3 



Equation (4) allows for an easy direct numerical access to the propagator once 
enough eigenfunctions, e.g., in position representation, are known. 



3 Wigner eigenfunctions and propagator for the quantum baker 
map 



The baker transformation [7] is arguably the simplest example of a fully chaotic 
map of the unit square to itself. A quantum version has been suggested in 
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Fig. 1. Wigner functions according to Eq. (1) for eigenstates of the quantum baker 
map, Eq. (5). Panels a,b show diagonal Wigner eigenfunctions W a ^ a for (a) a = 88 
and (b) 159. Panels c,d are real and imaginary parts, respectively, of the off-diagonal 
Wigner eigenfunction Ws8,i29- Hilbert-space dimension is N = 200. Greyscale ranges 
from negative (dark) to positive (bright) values. Unstable periodic orbits of the 
corresponding classical map underlying the scars have been marked with crosses. 

Ref. [7] in terms of Fourier transformations F M , (F M ) mjn = M~ l l 2 e~ 2mmn / M 
of sizes M = N and N/2, where N (even) denotes the Hilbert-space dimension, 
as follows, 



B 



N 



Fn/2 



F, 



(5) 



N/2 



so that wavefunctions map i/)'(n') = Z) n =o 1 (^A r )n'nV ; ( n )? an d the density op- 
erator accordingly. We have calculated diagonal Wigner eigenstates W a>a and 
off-diagonal ones W a> p and used them to construct the Wigner propagator 
equivalent to the unitary map (5) according to Eq. (4). The eigenstates a = 88 
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Fig. 2. Propagator of the Wigner function according to Eq. (4) for initial condition 
m' = (6, 25) (bold cross) on the period-3 unstable periodic orbit as marked in 
Fig. 1, after (a) 1 and (b) 2 applications of the same quantum baker map as in 
Fig. 1. Hilbert-space dimension is N = 42. Other parameters and greyscale as in 
Fig. 1. 

(Fig. la) and 129 (lb) clearly show scars of classical unstable periodic or- 
bits (crosses) which up to now could be identified only in the corresponding 
Husimi functions [9], with considerably inferior resolution. Figures 1 c,d show, 
respectively, the real and imaginary part of the off-diagonal eigenstate 1^88,129, 
scarred simultaneously by both orbits that appear in the eigenstates W^ss 
and Wi29,i29- 

The Wigner propagator is depicted in Fig. 2, for fixed initial m' = (6, 25) 
at t' = 0, the leftmost point of one of the period-3 orbits scarring ^129,129 
(Fig. lb), at final times t" = 1 (panel a) and t" = 2 (b). The strong localization 
of the propagator on the underlying classical orbit is evident, as are quantum 
coherence effects reflected in oscillatory patterns far off the classical orbit. 



4 Conclusion 

The redundance-free Wigner function for discrete phase spaces proposed here 
enables to identify classical and other phase-space structures in states and time 
evolution of quantum systems with finite-dimensional Hilbert space, without 
being marred by artificial ghost images. It facilitates developping and assessing 
semiclassical approximations for the propagator in such systems [10], provid- 
ing an alternative point of view for the analysis of coherent structures like 
"quantum carpets" [11,12]. 

The method allows to represent signals simultaneously in time and in fre- 
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quency space, in all applications where data are both discrete and periodic in 
time (or space) such as in optics, acoustics, or other fields possibly far away 
from quantum mechanics. 
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